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Abstract 

We study the stability of the ghost dark energy model versus perturbation. Since this kind of 
dark energy is instable in Einsteinian general relativity theory, then we study a new type of 
Brans-Dicke theory which has non-minimal coupling with matter which is called chameleon 
Brans-Dicke (CBD) model of gravity. Due to this coupling the equation of conservation 
energy is modified. For considering the stability of the model we use the adiabatic squared 
^ sound speed, Cg, whose sign of it determines the stability of the model in which for > 

Qi the model is stable and for < the model is instable. However, we study the interacting 

^ \ and non-interacting version of chameleon Brans-Dicke ghost dark energy (CBDGDE) with 

cold dark matter in non fiat FLRW metric. We show that in all cases of investigation the 
model is stable with a suitable choice of parameters. 
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1 Introductions 



The cosmological and astrophysical observations such as type la supernovae data[T], Wilkin- 
son Microwave Anisotropic Probe (WMAP) [2j, X-ray [3j, large scale structure ^ and ete, 
indicate that our universe is in accelerating expansion phase. People have introduced an 
energy component of the universe, called dark energy to describe this acceleration. The sim- 
plest model of DE is a tiny positive time-independent cosmological constant, A, for which 
has the equation of state a; = — 1 [21 E]- Plenty of other DE models have also been proposed 
for explain the acceleration expansion either by introducing new degree(s) of freedom or by 
modifying gravity [HElElIinillllllSlIiaiiaE]. 

In recent years there has been a new attention to the so called scalar-tensor gravity. 
The scalar-tensor models include a scalar field, cj), whit non-minimal coupling to the ge- 
ometry in the gravitational action, which has been introduced by Brans and Dicke (BD) 
|16j . They proposed a scalar degree of freedom to incorporate the Mach's principle into 
general relativity. The mechanism that creates a non-minimal scalar field coupling to the 
geometry, can also lead to a coupling between the scalar and matter field. Therefore authors 
|17j introduced a scalar field which it has a coupling to matter with order unity strength, 
named chameleon mechanism. Indeed, the chameleon proposal provides a way to generating 
an effective mass for a light scalar field via the field self interaction, and the interaction 
between matter and scalar field. 

Recently, a new kind of dark energy model, so called "ghost dark energy" has been 
investigated [181 119j. Originally, the Veneziano ghost was introduced as a solution to U{1) 
problem in low energy effective theory of QDF [20] . The ghosts make a small energy density 
contribution to the vacuum energy due to the off-set of the cancelation of their contribution 
in curved space or time-dependent background. The authors of [19] have clarified the 
decoupling of the QCD vector ghost to the vacuum energy density in the Rindler space- 
time. They have fund that it gives the vacuum energy density proportional to Hubble 
parameter, HAq^j^ of the right magnitude ~ (10~^et;)^, where H is the Hubble parameter 
and Aqcd'is QCD mass scale. The authors of [19] have climbed that in the ghost model of 
dark energy, one needs not to introduce any new degree of freedom or modify gravity and 
it is totally embedded in standard model of gravity. But according to results of |2U the 
GDE model in the Einsteinian theory of gravity has a behavior like a cosmological constant 
at the late time (w = —1) and the equation of state can never cross —1, and this is similar 
to the behavior of quintessence. Also, they have shown that the adiabatic squared sound 
speed of the GDE model in context of standard model of cosmology is negative and then 
this model can not be stable. Therefore, in this work we want to consider the stability of 
GDE model in context of BD model which has a non-minimal coupling with matter. The 
fundamental key quantity for studying the stability of a model is the squared adiabatic 
sound speed which is obtained a small perturbation in the back ground energy density [22j. 
The sign of plays a crucial role in determining the stability of the background evolution. 
If cl < 0, it means that the model is classically instable against perturbation. This issue 
has already been investigated for some DE models such as chaplygin gas and tachyon DE 
|23j , holographic DE [24| , agegraphic model of DE [25j . 

In this work, we study the ghost model of dark energy in the context of chameleon 
Brans-Dicke model of gravity. We investigate the cosmological evolution of our model 
with/without interaction between DE and cold dark matter. We analytically and numer- 
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ically compute some quantities such as scale factor a, EoS parameter of dark energy uJd, 
deceleration parameter q, fraction of dark energy ^l^, squared adiabatic speed of sound 
and so on. 

This work is organized in four sections, of which this introduction is the first. In section 
two, the action is introduced and the field equation, the scalar field equation of motion, 
the modified conservation of density energy are obtained. In section 3 the interacting and 
non interacting GDE model of CBD theory is investigated in non flat FRW space-time, and 
section 4 is the summarize of our results. 



2 General Framework 

For our investigation, we consider the chameleon-Brans-Dicke action 



S = J d'x^ (^R - '^d^cPdf'cP + 2fi(P)Cm^ 



where g is the metric determinant, R is the Ricci scalar constructed from the metric 5^,^, and 
(l)R has been replaced with the Einstein-Hilbert term is such a way that G~jjr = 167r(p, cj) is 
the chameleon-Brans-Dicke scalar field, oj is the dimenssionless Brans-Dicke constant. The 
last term on the right hand side of ([T|), f{(l))Cm, indicates non-minimal coupling between the 
scalar field and matter. One can obtain the gravitational field equation by taking variation 
of the action ([1]) with respect to the metric g/^n, 



R 



1 J^_Sy^rj. ,rj.4> 



(2) 



where Tn,^ indicates the scalar field energy-momentum tensor which is 



rp<f> 



OJ 



1 



+ 



1 



(3) 



here □ is the four dimensional d'Alambert operator, and T^^ indicates the matteiH energy- 
momentum tensor which is defined 



2 5{^Cr 



and represented by 



6g. 



(4) 



(5) 



Where Un is the four-vector velocity of the fluid satisfying UnU^ 



-1, pt and pt are 



respectively the total energy density and total isotropic pressure of the barotropic perfect 
fluids which have fllled the universe. Taking variation of action with respect to scalar fleld 
(p gives us the Klein- Gordon equation for the scalar fleld as 



{2u + 3)n<p = [/(0)r - 2<pf'{<p)Cr 



(6) 



^ In fact matter stress-energy tensor consists of all perfect fluids stress-energy tensor, namely T^^ — 
T^J -\- T^uf^^ + T^l) + T^i-' . Here the subscript b, cdm, r and d indicate baryonic matter, cold dark matter, 
radiation and dark energy respectively. Indeed in this work, the same as others, we assume that dark energy 
has an averaged bahavior like perfect fluid. 
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where, Lm is the Lagrangian of the matter and T is the trace of matter stress-energy tensor. 
The Bianchi identities, together with the identity (DVa — Van)!^ = ^ah^^'^c-, imply the 
non-(covariant) conservation law 

VaT»^= [5'^''£„ + r»'']^, (7) 

and, as expected, in the limit /(</>) = constant, one recovers the conservation law SI J""^ = 0. 

Our aim in this work is to consider the ghost model of dark energy in context of CBD 
model in the Friedmann-Robertson- Walker (FRW) Universe which is described by the fol- 
lowing line element 

ds^ = dt' - a{tf ( + r'dnA , (8) 
\ 1 — kr^ J 

where a{t) is the scale factor, and k is the curvature parameter with k = —1,0,1 corre- 
sponding to open, flat, and closed Universes, respectively. By making use ([2]), ([5]), ^ 
and dHl) one can arrive at 



2H + 3H = - p,- 2H{-)--, (10) 



4> + ^H4, = (2a; + 3) [fi4>)T-2ct>f'i^)Cm], (11) 

where H = d/a is the Hubble parameter. In the above equations, the EoS parameter of the 
baryonic and dark matter is pt = uJtPt- 

In |26j . it is shown that a natural choice for the matter Lagrangian density for perfect 
fluids which based on ([3]) can give us the stress-energy tensor, ([5]), is Cm = Pt, where pt is 
the pressure. However, although Cm = Pt does indeed reproduce the perfect fluid equation 
of state, it is not unique. For example, other choices are Cm = —pt or —ntUt, where pt is 
the energy density, nt is the total particle number density, and at is the total physical free 
energy defined as at = pt/nt — Ts, with T being the fluid temperature and s the entropy 
per particle [26]. Therefore one may introduce the matter Lagrangian a^ 

1 3 

= -^Pt + -^Pt, (12) 

which can give us the stress-energy tensor, based on ([3D. So using ([Sj) we can rewrite 
the tt component of ([7]) as follows 

pt + 3H{l + ujt)pt = -^jil+uJt)pt. (13) 

It is seen that there is an addition parameter in evolution equation of model. So ac- 
cording [27j we shall assume that the scalar field can be introduced as a power law of the 
scale factor as 

^ = Ar(^)« = iVa«, (14) 



^As mentioned before, the choice of Cm is not unique but only for simphcity we have chosen it as (|12p 
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here ao is the scale factor at the present time, ao = 1, and N = constant. In fact there is 
no compehing reason for this choice. However, it has been shown that for small it leads 
to consistent results and the product |^|a; results of order unity [27]. According to the solar 
system experiments the magnitude of oj is more than 40000 (w > 40000) [28' '29] and this 
means that the value of ^ has to very small ~ lO"''. Using the latest WMAP and SDSS data, 
the observational constraints on BD model in a flat Universe with cosmological constant 
and cold DM is obtained [3U]. They found that within 2a range, the value of u satisfies 
(jj < —120.0 or OJ > 97.8. They also obtained the constraint on the rate of change of G at 
present as -5.53 x 10"^°s"^ < G/G < 3.32 x 10"^°s~\ at 2a confidence level. So in our 
case with assumption 1^ we get -5.53 x lO-^o^-i < g/g = ^/(j) = < 3.32 x lO-^Og-i. 
This relation can be used to put an upper and lower bound on ^. Assuming the present 
value of the Hubble parameter to be Hq ~ 2.11 x 10~^^s~^, we obtain 

- 0.026 < C < 0.015. (15) 

For getting a better insight, we continue our work based on a power law form for f{(j)) 
as /(0) = /o<A*, where /o and here s = are constant. Also we assume there are only 
two components GDE and CDM in the Univers^. Therefore pt = Pm + Pd = Pd and then 
Pt = oJtPt = ^dPd, where pt = Pd + Pm- Hence from (10) and using (11), we have 

p„, + 3Hil + ^)pm = 0, (16) 
Pd + 3i/(l+Wd)(l + |)pd = 0. (17) 

One can interpret the interaction term in (|13p and also in (jl6p as pm = (,pm/^- This 
means that for ^ < 0, the interaction between chameleon scalar filed and matter fluid can 
create a negative pressure. Therefore we expect this process help to describe the positive 
accelerating expansion of Universe. 



3 GDE in the context of CBD model of gravity 

In this section we consider the GDE in a non flat space-time of FRW Universe in chameleon 
BD model of gravity. The ghost energy density is proportional to the Hubble parameter 

m- 

Pd = aH, (18) 
where a is a constant. According to the results of [19j, a ~ 

^QCD' where A ~ lO^MeV 

is QCD mass scale. Therefore the value of dark energy at the present time, with H ~ 
10-3^MeV is about (3 x IQ-^eV)^. This value is in an excellent agreement with observed 
DE density [19] . As mentioned in the Introduction this model of dark energy in the context 
of standard model of cosmology is not stable and this means that GDE model in standard 
model of gravity has shortcoming and need to study in another model of gravity. 

^The results of this study based on s 7^ 1 is not very different to tfie our investigation. 

^ Since we are interesting to study tlie our model at the late time and in this stage dark energy is 
dominant, then for simplicity we absorb other components of perfect fluid(baryonic matter and radiation) 
in cold dark matter part. 
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3.1 Non interacting GDE 

In this subsection we investigate the GDE model in the CBD framework in a non flat FRW 
space time and we will obtain the EoS parameter, deceleration parameter, the evaluation of 
fractional energy density and the adiabatic squared sound speed of GDE. Taking the time 
derivative of equation Q, using relation (jl4p and the continuity equations (|17p . we find 

= -3/3(1 



.l|3/30_(2-3/3)J7fc|- A^A, (19) 



where 6 = 1 — ^(a;^/6 — 1), /? = l + s^/4. One can obtain the EoS parameter of GDE model 
of chameleon BD as 

^ ^ ^A: ^20) 

28 — r2(^ 26 — ilj; 

where = (2 — 3/3)0^/3/3. It is also interesting to study the behavior of the deceleration 
parameter defined as 

. H 

{3/3-l)nd ^ (3/3-2)(g + 17fc) ^2^^ 
26 — 26 — 

Note that ^ is very small and ^ S (—0.026,0.015). Hence according to definition of 6, it 
must be in (—3.53, 1). This means that 6 can accept the negative value. But, observational 
data indicates that the current expansion of Universe is accelerating. Several attempts have 
been made to justify this accelerated expansion and this constraint requires ujd < —1/3 and 
q < 0. Therefore, from (j2ip one can see that, this model can explain the positive accelerating 
expansion of Universe if 6 > $1^/2- Then the physical value of 6 is in ^Idf^ < 6 < 1 interval. 
Also, the equation of motion of GDE can be obtained by using (flSj) and (JT9]) as 

^, 2^_ai ^ 3^^^ ^ j^^j _ 3^^^^ ^22) 

here prime denote derivative with respect to efolding number x = In a. Integrating of ([22] ) 
gives 

26ln{nd) - AlnilkQd - 1|) = 3pAx + c, (23) 

where 

A = 6-nk 
1 



k 



and 

c = 26ln{^}do)-Aln{\kndo - M) 
where is the present value of dimensionless energy density of DE. 

At this stage we consider the stability of GDE in the context of CBD model of gravity. 
The main key quantity for investigating the stability of a model is the squared adiabatic 
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sound speed which is obtained a smah perturbation in the back ground energy density. To 
derive the relevant synchronous-gauge equations of motion, the variables that characterize 
the fluid are linearized about a spatially homogeneous background |22] 



p{t,x) = pb{t) + 5p{t,x), 
p{t,x) = pb{t) + c'^Jp{t,x). 



(24) 
(25) 



The subscript b denotes the spatially homogeneous background (mean) value of the corre- 
sponding quantity, and 

is the squared adiabatic sound speed of the fluid. Using energy conservation equation yields 

m 

1 92 



V'6p{t,x)--—^[6p{t,x)]=0. 
c„ 01 



(27) 



This equation is an ordinary wave equation. When > the answer of ()27p is as Sp{t, x) = 
5po exp(±i[a;t-|-A;x]), which is an oscillatory waves and this indicates a propagation mode for 
the density perturbations. But, when < 0, in this case the oscillations becomes hyperbolic 
and the density perturbations will grow with time as 6p{t,x) = 5pQeKp{ujt ± ikx). Thus 
the perturbation is growing and this means the background energy of the model is instable. 
Therefore, as was mentioned before, the main key quantity for investigating the stability of 
a model is the squared adiabatic sound speed. So by using (p6]) and EoS pd = (^dPd we have 



2 P , . Pd 

P Pd 



taking time derivative of ()20p . we have 



and then we have 



(^d— = {0 - ^k) 
Pd 



{29 - ^d? 

nd-e 



(29 - ^d? 



(28) 



(29) 



(30) 



since VL^ < 0, then it is clearly seen that for 9 < Cld the squared adiabatic sound speed is 
positive and then the model is stable. 



3.2 Interacting GDE 

Some observational data such as, observational of the galaxy cluster Abell A586, supports 
the interaction between DE and DM |25) . Therefore in this section we introduce the direct 
interaction between DE and CDM and study the evolution dynamics of the model. Therefore 
according to (fT6|) and (fT7|l the conservation equations are modified as 

p^ + mpp^ = Q, (31) 
Pd + SHpd(3{l+Ud) = -Q, (32) 

where Q denotes the interaction between GDE and CDM. A generic form of Q is not avail- 
able. Three forms which are often discussed in the literature are asQ = Sb^Hpd, 3b^Hpm, Sb^Hpt, 
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where 6^ > 0. We want investigate the model for these three kinds of interaction. So that 
we consider it with a more general procedure with respect to the pervious subsection. Hence 
we rewrite ([9]), ([31]) and dMI) as 

Qm + ^d = e + Qk, (33) 

nm + —nm + 3(3Hnm = (34) 

(led+ ^^ed + m'^+^d)H^d = (35) 

where = Pm/^H^, ^d = a/3H, and il^ = k/SH'^a^. Using (30) and (31) we have 

+ —{9 - ^d) + 3PH{e -nd) = ^ + {2- ^m^k, m 



and by putting (j35|) in (j36j) one can get 

2eH + 3/3uJdH^nd + 3/3H^e = (2 - 3/3)^/^^^. (37) 
Note that in GDE model we have 

HQd = HqQ,o, (38) 

where Hq and Qq are the value of Hubble parameter and fraction of the ghost dark energy 
at the present time. Expressing (I36p in terms of efolding-number a; = In a, and making use 
(1381) we have 

" ^^^^^ = MP~ ^^^^ " ^'^^ + " ^^^^'^^ ^^^^ 
Here, we can obtain the equation of state parameters of the GDE versus VL^ as 

'"'^ = 3/3(20 - 0.) ' 

where = Q/{3Q.dH^) and deceleration parameter is as 

(Oq + 3/3 - l)^lrf (3/3-2)(r?fc + g) 
^ (20 -f^,) + {26 -Qd) ' ^ ^ 

Although (3/3-2) > 0, but since (3/3-2)(6l + ilfc) is smaller than (rig + 3/3 - 1)17^, therefore 
the deceleration parameter of the GDE model in BD scenario is negative for 6 > Vtd/2, 
and one can find that the EoS parameter of GDE in our model, w^, can cross the phantom 
divide line, ujd = for 

17, 3P{^d + ^k) . 

T<^< 3/3-2f]Q • 

Also we can obtain the squared adiabatic sound speed of the our model as 

2 _ dpd 



dpd 



l-Qd 



ojd. (43) 



For getting better insight we consider the interacting GDE for three different forms of Q 
below. However for the sake of briefness, we will investigate the case of Q = Sb^Hpd in 
detail. 
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Figure 1: (a): This sub-figure shows versus for b = 0.25, ilk = 0.02 and different 
values of C= (0.0063, red(solid)), (0.0068, black(dashed)), (0.0073, blue(dotted)), (0.0085, 
green (dashed-dotted)). (b): This sub-figure shows versus ild for 9 ~ 0.63, = 0.02 
and different values of 6=(0, red (solid)), (0.1, black(dashed)),(0.25, blue (dotted)), (0.45, 
green (dashed-dotted)). We have taken Q = 3b^Hpd. 




0.000 0.002 0.004 0.006 0.008 
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0.000 0.002 0.004 0.006 0.008 

(b) 



Figure 2: (a): This sub-figure shows Ud versus ,f for 6 = 0.25, Qk = 0.02 and Qd = 0.76. 
(b):This sub-figure shows q versus for b = 0.25, ^7^0. 02 and Qd = 0.76. We have taken 
Q = 3b^Hpd. 

3.2.1 Q = Sb'^Hpd 

In this case Qq = 3b^ and equations (00]) and (jH]) reduce to 

9 Clk 29b'^ 

~ ~ {29 - Qd) ~ {29 - Qd) ~ mo - ^d) ' ^ ^ 

(36^ -^ 3/3 - l)l]rf {2,p-2){ilk + 9) 
^ {29 -^d) ^ {29 -^d) ' ^ ' 



Although ilfc < but the magnitude of it is smaller than other terms, therefore (j44p shows 
that ojd is always negative. We plot the EoS parameter versus fi^ in figure 1. Figure 
l.a shows the behavior of a;^ versus Vid for b = 0.25, il^ = 0.02 and different values of ^= 
(0.0063, red(solid)), (0.0068, black(dashed)), (0.0073, blue(dotted)), (0.0085, green(dashed- 
dotted)). Note that 9 € (0.52,0.74) for this values of It is clearly seen that oJd cross the 



9 




Figure 3: (a): This sub-figureshows q versus 0,^ for b = 0.25, fijt = 0.02 and different 
values of ^= (0.0063, red(solid)), (0.0068, black(dashed)), (0.0073, blue(dotted)), (0.0085, 
green (dashed-dotted)). (b):This sub-figure shows Ud versus for 9 ~ 0.63, = 0.02 
and different values of 6= (0, red (solid)), (0.1, black(dashed)),(0.25, blue(dotted)), (0.45, 
green (dashed-dotted)). We have taken Q = Sb'^Hpii- 



line ojfi = —1/3 at nearly Q^, = 0.35 and also cross the phantom divide line (w^ = —1) at 
almost Qd G (0.58,0.77) in which the crossing point is different via the magnitude of ^. 
Figure l.b indicates Ud versus i^d for 6 ~ 0.63(^ = 0.0073), fi^ = 0.02 and different values 
of b. It is seen that for 6 = (red-solid) the evolution of Universe inter to the positive 
accelerating expansion at nearly 0^ ~ 0.36 and cross the phantom line at ~ 0.64. This 
figure shows that due to the direct interaction between DE and matter the phantom line 
crossing take place earlier. Also we plot the EoS parameter and deceleration parameter of 
GDE versus ^ in figure2. This figure shows that for all ^ G (0, 0.0094) the evolution of the 
Universe completely is in positive accelerating expansion phase and uid can cross the line 
ujd = -1- 

Also we plot the deceleration parameter, q, versus ild in Fihure3. The figure 3. a shows 
the behavior of deceleration parameter for b = 0.25, = 0.02 and different values of ^. 
According to this figure, the evolution of Universe inter to the accelerating expansion phase 
at nearly ild — 0.35 for all value of ^ G (0.0025,0.0094). Figure3b, shows q versus il^ for 
9 ~ 0.63 (^ = 0.0073), Qk = 0.02 and different values of b. This figure indicates that the 
evolution of Universe inter to the positive accelerating expansion for — 0.3 for large 
interaction between GDE and matter which is not agree with the stability of the model. 
But the evolution process of Universe inter to the accelerating phase at nearly ^Id — 0.37 
for 6 = (without interaction) and this is in a good agreement with stability of model, 
because the stability is satisfied only when ild/2 < 9 < i^d- 
In this case the adiabatic squared sound speed is 



+ Ofc + — 6^ 



(46) 



(20 - SU) 



In this equation 9, fS and also 6^ are positive and although (ij^ is negative, but the value 
of it is very smaller than the other terms, so the squared sound speed is positive only for 
r2(i/2 < 9 < 0,d- This means that this interacting version of model can be stable. We plot 
the adiabatic squared sound speed of GDE versus ^Id in figure 4. This figure explicitly show 
that for some values of fi^, especially at the present, the squared sound speed of GDE of 
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Figure 4: (a): This sub-figure shows versus Qd for b = 0.25, = 0.02 and different 
values of ^= (0.0063, red(sohd)), (0.0068, black(dashed)), (0.0073, blue(dotted)), (0.0085, 
green (dashed-dotted)). (b): This sub-figure shows uj^ versus ild for 9 ~ 0.63, = 0.02 
and different values of 6= (0, red (solid)), (0.1, black(dashed)),(0.25, blue(dotted)), (0.45, 
green (dashed-dotted)). We have taken Q = Sb'^Hpd- 



CBD model of gravity is positive. Figure 4a indicate against 0,^. for b = 0.25, fl^ = 0.02 
and different values of ^ and shows that can be positive for Q^. ^ 0.55. Figure 4b presents 
versus ^l^ for 6 ~ 0.63(^ = 0.0073), il^ = 0.02 and different values of b for the first forme 
of interaction, Q = Sb^Hpd- This figure tells us that a direct interaction between GDE and 
matter, such as Q = Sb^Hpd, in the cotext of CBD model of gravity has not any effect on 
the stability of the model. Because for all value of b, even 6 = (without interaction), the 
adiabatic squared sound speed has the similar behavior versus Orf and > at Qd ~ 0.63. 
This means that the GDE of CBD model is stable at the present time. The relation Cg for 
an interacting case of GDE with matter is plotted versus ^ in figure 5. This figure shows 
that for interval ^ G (0.0061,0.0094) which is equivalent to 6^ G (0.42,0.745) the adiabatic 
squared sound speed of the model is positive and then the model is stable in this form of 
interaction. 

Figure 5 present versus Cg is plotted for Q = 3b^Hpd in figure 5a and for all three 
forms of Q in figure 5b. These two sub-figures show that the behavior of versus ^ is 
similar for all three forms of interactions. 
Finally the equation (j39p becomes 

^,^2^-J^, ^ ^^^^ _ ~ _ ^^^2 ^ ^^^^ 

its analytical solution reads 

201nl7d-Oiln(|l- fciOrfl) = 3/3(0 - fifc)x + c, (48) 

where 

«. ^ ,49, 
= ^±^. (50, 

13(0 -nt) 
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Figure 5: (a):This sub-figure shows versus ^ for 6 = 0.25, ilfc = 0.02 and 17^ = 0.76 in 
the case of Q = Sb'^Hpii. (b): This sub-figure shows versus ^ for 6 = 0.25, 17^0.02 and 
Qd = 0.76 in the case of Q = 36'^ifpd(yellow), Q = 36^ -ffpm (green), Q = 36^ ffpt (black). 



and c = 201n$7rfo ~ fiilndl — kiQ^ol) is the integration constant and is the fraction 
of dark energy at present time. The relation of versus x = ln(a) is shown in Figure 6. 
From the figure6 we can see that varies from at early time to 1 at late time. 




X 



Figure 6: This figure show Qd versus x = ln(a) for 6 = 0.1, ujk = 0.02 and different values 
of C = (0.0065, black(dashed)), (0.0075, green(dotted)), (0.0085, red(solid)). We have taken 



3.2.2 Q = Sb^Hpm 

In this case Qq = 3b'^{6 + ftk — ^d)/^d and equation (l40l) . (j4T]) and (jiSl) reduce to 

je + nk) 2eb^ 2eb''{e + ^k) . . 

'^'^ {26 - nd) me - ^d) P^d{2e -^dY ^ ^ 

(3/j _ 3^2 _ 1)^^ (3/3-35^-2)(g + Ofc) 

^ {26 -nd) ^ {26 -nd) ' ^ ' 

- ^^\2e-ndY^'''nd{26-ndr ^ ^ 
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fi_d D._d 

(a) (b) 

Figure 7: (a): This sub-figure shows versus Jl^ for b = 0.25, 0^ = 0.02 and different 
values of = 0.0063 (red-solid), 0.0068 (black-dashed), 0.0073 (blue-dotted) in the case of 
Q = 3b^Hpm- (b): This sub-figure shows ojd versus Cld for b = 0.25, = 0.02 and different 
values of ^ = 0.0063 (red-solid), 0.0068(black-dashed), 0.0073( blue-dotted) in the case of 
Q = Sb'^Hpt. 



where 

29 

uji = 9 + nk- jb^, 

W2 = ^eb\e + nk). 

Prom the above equations are seen that when — 0, ujd, Qi tend to —b'^{6 + Q,k)/P^di 
(3/3 — b^ — 2){9 + r2fc)/20, —uj2/9^d respectively. We see that at the early time and Cg 
are divergent. 

The equation of motion for il^ is 



(P-b^)e-(--i3 + b')ak\ - 3(13 - b^)aj. 



Its analytical solution gives 

29lnnd - n2ln{\l - ki^dl) = \3{f3 - b'^)9 - {3b'^ + 2 - 3/3)^1^ 

where 

(/3 - b^)9 + {b^-(3 + 2/3)VLk 



x + c, 



k2 



P-b"^ 



{l3-b^)9-{b^- /3 + 2/3)nk' 
and c = 29lni^fiQ — r22ln(|l — k2^do\) is the integration constant. 



(54) 

(55) 

(56) 
(57) 
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0.0 0.2 0.4 0.6 0.8 1.0 

n_d 

(a) 



0.0 0.2 0.4 0.6 0.8 1.0 

n_d 

(b) 



Figure 8: (a): This sub-figure shows versus Ud for b = 0.25, Qk = 0.02 and different 
values of C = 0.0063 (red-soUd), 0.0068 (black-dashed), 0.0073 (blue-dotted) in the case of 
Q = 3b^Hpm- (b): This sub-figure shows versus 17^ for b = 0.25, 17^ = 0.02 and different 
values of ^ = 0.0063 (red-solid), 0.0068(black-dashed), 0.0073( blue-dotted) in the case of 
Q = Sb'^Hpt. 

3.2.3 Q = Sb^Hpt 

In this case nq = 36^(61 + nk)/^d and eguationdiOl). (l4T]l and dS]) reduce to 

(9 + nk) 2ebH9 + ^k) 



(29 -n^) (3nd{29-ndy 

(3/3-l)J^d (3/3-36^ -2) (e + f^A 



2u> 



{29 - Qd) 
Qd-O 



+ 



+ UJ2 



{29 - nd) 

SQd - 4(9 



where 



{29-ndy 'nd{29-Qd?'' 



uj\ — 9 -\- r^jij. 



(58) 
(59) 
(60) 



U2 = ^9b\9 + nu). 

It is seen that the quantity Wd, g, have the same behavior which we obtain for Q 
Sb'^Hpm- This means that cod and are divergent in the early time. 



Finally in this case the equation of motion for dimensionless density energy parameter 



IS 



1 29 — 



{p-b^)9-{--/3 + b^)nk -3i3nd- 



Integration (l47l) . gives 

20 In fid - ^3 ln(|l - feafia 

where 



3(/3 - b'^)9 - (36^ + 2 - 3/3)fiA 



k3 



{f3 - b^)9 + {b^-f3 + 2/3)fifc 

^ 

{/3-b^)9-{b^- f3 + 2/3)nk' 



(61) 

X + C3, (62) 

(63) 
(64) 
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and 

C3 = 201nOdo-5^3 ln(|l-fc3f^dol)- 

We plot ujd versus Qd for Q = 3b^Hp„i and Q = Sb^Hpt in figure 7. These two sub-figures 
show that the behavior of uj^ versus 0^ is similar for two different kinds of interaction and 
also they show the evolution of Universe is completely in positive accelerating expansion 
phase, because Ud is always less than —1/3. These figures indicate ujd will increase from 
— oo to a local maximum below ujd = —1/3 and then decrease to the below of Ud = — 1- 
This means that Ud cross the line = — 1 two times. One time from phantom phase 
to quintessence phase nearly at early time and another time from quintessence phase to 
phantom phase in which at the present {^Id = 0.76) the second crossing has took place. 

Figures 8(a) and 8(b) show the squared sound speed versus Qd for Q = 'ib^Hp^ and 
Q = 3b^Hpt respectively. These figures indicate is negative in < $7^ < 0.63 for these 
two kinds of interaction, and for > 0.63 one can obtain a positive value for by a 
suitable choice of ^ parameter. This means that for these two kinds of interaction the GDE 
model is stable in the present time in the context of CBD model of gravity. 

4 Conclusion 

Recently, Ghost dark energy is introduced [19]. This kind of dark energy is a phenomeno- 
logical vacuum energy which is rooted from the Veneziano ghost of QCD. The ghost dark 
energy model is proportional to Hubble parameter. According to results of [21j, the adia- 
batic squared sound speed of the GDE model in the context of Einsteinian theory of gravity 
is negative and then this model is instable. And also it is shown that this model has a 
behavior like a cosmological constant at the late time (w = —1) and the equation of state 
can never cross —1. We study the stability and the evolution of this model in the context of 
Brans-Dicke model which has non-minimal coupling with matter, namely chameleon Brans- 
Dicke (GDECBD) model. At first, we investigated the GDE of CBD model without any 
interaction between GDE and CDM in a non flat FLRW metric. In this investigation we 
obtain the EoS parameter, deceleration parameter, the equation of motion for GDE fraction 
parameter and the adiabatic squared sound speed of the model. The obtained quantities 
show that the GDE model in the context of CBD scenario can describe the positive accel- 
erating expansion phase of Universe and the EoS parameter of dark energy can cross the 
phantom divide line, also the squared sound speed of the DE is positive, for a suitable choice 
if ^ parameter. We studied the cosmological dynamics of the model by considering three 
usual forms of interaction between DE and CDM. Considering the evolution of dimension- 
less density energy parameter, we have shown that Qd varies from at early time to 1 at late 
time. We fund that the evolution of Universe in GDE model of CBD scenario is completely 
in positive accelerating expansion phase for three forms of interaction. Also we obtained 
LOd crosses the line Ud = from phantom phase to quintessence phase for Q = Sb'^Hpd 
case, but the behavior of LOd versus Q,d is similar for Q = Sb^Hpm and Q = 3b^Hpt and also 
ojd cross the line cj^ = — 1 two times. One time from phantom phase to quintessence phase 
nearly at early time and another time from quintessence phase to phantom phase in which 
at the present {Qd = 0.76) the second crossing has took place. 

Finally we found the squared sound speed of GDE model in CBD framework for three 
kinds of interaction in the FRW Universe. All of our calculations show that is negative 
in < rjrf < 0.63 and for Qd > 0.63 one can obtain a positive value for by a suitable 
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choice of ^ parameter. This means that the interacting and non-interacting case of GDE 
model is stable in the present time in the context of CBD model of gravity. 
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